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Abstract Using the recent formulation of Noether's theo- 
rem for the problems of the calculus of variations with frac- 
tional derivatives, the Lagrange multiplier technique, and the 
fractional Euler-Lagrange equations, we prove a Noether- 
like theorem to the more general context of the fractional op- 
timal control. As a corollary, it follows that in the fractional 
case the autonomous Hamiltonian does not define anymore 
a conservation law. Instead, it is proved that the fractional 
conservation law adds to the Hamiltonian a new term which 
depends on the fractional-order of differentiation, the gener- 
alized momentum, and the fractional derivative of the state 
variable. 
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vation laws in the calculus of variations and optimal control 
is to reduce the number of degrees of freedom, and thus re- 
ducing the problems to a lower dimension, facilitating the 
integration of the differential equations given by the neces- 
sary optimality conditions. 

Emmy Noether was the first who proved, in 1918, that 
the notions of symmetry and conservation law are connected: 
when a system exhibits a symmetry, then a conservation law 
can be obtained. One of the most important and well known 
illustration of this deep and rich relation, is given by the 
conservation of energy in Mechanics: the autonomous La- 
grangian L(q,q), correspondent to a mechanical system of 
conservative points, is invariant under time-translations (time- 
homogeneity symmetry), ancQ 



d 
dt 



[L(q,q)-d 2 L{q,q)-q] =0 



(1) 



o 

£>■ ■ 1 Introduction 

V , 

. The concept of symmetry plays an important role both in 
C3 Physics and Mathematics. Symmetries are described by trans- 
formations of the system, which result in the same object 
after the transformation is carried out. They are described 
mathematically by parameter groups of transformations. Their 
importance ranges from fundamental and theoretical aspects 
to concrete applications, having profound implications in the 
dynamical behavior of the systems, and in their basic quali- 
tative properties. 

Another fundamental notion in Physics and Mathematics 
is the one of conservation law. Typical application of conser- 
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follows from Noether's theorem, i.e., the total energy of a 
conservative closed system always remain constant in time, 
"it cannot be created or destroyed, but only transferred from 
one form into another". Expression (03 is valid along all the 
Euler-Lagrange extremals q(-) of an autonomous problem of 
the calculus of variations. The conservation law ([TJ is known 
in the calculus of variations as the 2nd Erdmann necessary 
condition; in concrete applications, it gains different inter- 
pretations: conservation of energy in Mechanics; income- 
wealth law in Economics; first law of Thermodynamics; etc. 
The literature on Noether's theorem is vast, and many ex- 
tensions of the classical results of Emmy Noether are now 
available for the more general setting of optimal control (see 
Il7l l29ll30l and references therein). Here we remark that in 
all those results conservation laws always refer to problems 
with integer derivatives. 

Nowadays fractional differentiation plays an important 
role in various fields: physics (classic and quantum mechan- 
ics, thermodynamics, etc), chemistry, biology, economics, 
engineering, signal and image processing, and control the- 
ory i2l ll3l[T6l . Its origin goes back three centuries, when in 
1695 L'Hopital and Leibniz exchanged some letters about 

1 Following 1111 . we use the notation <9,L to denote the partial 
derivative of function L with respect to its z'-th argument. 
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the mathematical meaning of -j± for n = 5. After that, many 
famous mathematicians, like J. Fourier, N. H. Abel, J. Liou- 
ville, B. Riemann, among others, contributed to the develop- 
ment of the Fractional Calculus ["13., 2011271. 

The study of fractional problems of the Calculus of Vari- 
ations and respective Euler-Lagrange type equations is a sub- 
ject of current strong research. F. Riewe [ 25 26] obtained a 
version of the Euler-Lagrange equations for problems of the 
Calculus of Variations with fractional derivatives, that com- 
bines the conservative and non-conservative cases. In 2002 
O. Agrawal proved a formulation for variational problems 
with right and left fractional derivatives in the Riemann- 
Liouville sense [1|. Then, these Euler-Lagrange equations 
were used by D. Baleanu and T. Avkar to investigate prob- 
lems with Lagrangians which are linear on the velocities 
[5|. In 11411151 fractional problems of the calculus of vari- 
ations with symmetric fractional derivatives are considered 
and correspondent Euler-Lagrange equations obtained, us- 
ing both Lagrangian and Hamiltonian formalisms. In all the 
above mentioned studies, Euler-Lagrange equations depend 
on left and right fractional derivatives, even when the prob- 
lem depend only on one type of them. In JT7 1 problems de- 
pending on symmetric derivatives are considered for which 
Euler-Lagrange equations include only the derivatives that 
appear in the formulation of the problem. In [8 9] Riemann- 
Liouville fractional integral functionals, depending on a pa- 
rameter a but not on fractional-order derivatives of order a, 
are introduced and respective fractional Euler-Lagrange type 
equations obtained. More recently, the authors have used the 
results of U) to generalize the classical Noether's theorem 
for the context of the Fractional Calculus of Variations II 1 1 11 - 
Differently from Hill , where the Lagrangian point of view 
is considered, here we adopt an Hamiltonian point of view. 
Fractional Hamiltonian dynamics is a very recent subject but 
the list of publications has become already a long one due to 
many applications in mechanics and physics [4 6 10. 12211231 
24,28]. We extend the previous optimal control Noether re- 
sults of [29 30 1 to the wider context of fractional optimal 
control (Theorem l4.2b . This is accomplished by means (i) of 
the fractional version of Noether's theorem [11], (ii) and the 
Lagrange multiplier rule J3J . As a consequence of our main 
result, it follows that the "total energy" (the autonomous 
Hamiltonian) of a fractional system is not conserved: a new 
expression appears (cf. Corollary |4.1b which also depends on 
the fractional-order of differentiation, the adjoint variable, 
and the fractional derivative of the state trajectory. 



2 Fractional Derivatives 

We briefly recall the definitions of right and left Riemann- 
Liouville fractional derivatives, as well as their main prop- 
erties im i20ll27l . 

Definition 2.1 Let / be a continuous and integrable func- 
tion in the interval [a,b]. For all f G [a,b], the left Riemann- 
Liouville fractional derivative a D"f(t), and the right Riemann 



Liouville fractional derivative ,D®f(t), of order a, are de- 
fined in the following way: 

» D ^" = r(^)(s)"/ ( '- ,,r °" /(9) '"'- ,2) 



o W'> = ra(-0T (9 -' r ""' /(e) '' ,, 



(3) 

where n 6 N, n— 1 < (X < n, and F is the Euler gamma 
function. 

Remark 2.1 If a is an integer, then from (O and (O one 
obtains the standard derivatives, that is, 

«/w= (-!)"/«. 

Theorem 2.1 Let f and g be two continuous functions on 
[a,b]. Then, for all t £ [a,b], the following properties hold: 

1. for p>0, 

a Df (/(f) +g(t)) = a D?f(t) +aD?g(t) ; 

2. for p > q > 0, 

M ( a D; q f{t))= a Dr q f{t); 

3. for p>0, 

X(«A~'7(0)=/(0 

(fundamental property of the Riemann-Liouville fractional 
derivatives). 

Remark 2.2 In general, the fractional derivative of a con- 
stant is not equal to zero. 

Remark 2.3 The fractional derivative of order p > of func- 
tion (t — a) v , D > 1, is given by 

a Df (t -a) v = F(U + 1) (t - a) v ~P . 

,v 1 r(- P +v + iy ' 

Remark 2.4 When one reads "Riemann-Liouville fractional 
derivative" in the literature, it is usually meant (implicitly) 
the left Riemann-Liouville fractional derivative. In Physics, 
t often denotes the time-variable, and the right Riemann- 
Liouville fractional derivative of f(t) is interpreted as a fu- 
ture state of the process /(f)- For this reason, right deriva- 
tives are usually neglected in applications: the present state 
of a process does not depend on the results of the future 
development. Following [3], and differently from 1111 . in 
this work we focus on problems with left Riemann-Liouville 
fractional derivatives only. This has the advantage of simpli- 
fying greatly the theory developed in 1111 . making possible 
the generalization of the results to the fractional optimal con- 
trol setting. 

We refer the interested reader in additional background 
■ on fractional theory, to the comprehensive book l27ll . 
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3 Preliminaries 

In IT] a formulation of the Euler-Lagrange equations is given 
for problems of the calculus of variations with fractional 
derivatives. 

Let us consider the following fractional problem of the 
calculus of variations: to find function q(-) that minimizes 
the integral functional 

/[<?(•)]= [ b L(t,q(t), a D?q(t))dt, (4) 

■la 

where the Lagrangian L : [a,b] x W x W — > R is a C 2 func- 
tion with respect to all its arguments, and < (X < 1 . 

Remark 3.1 In the case a = 1, problem Q is reduced to the 
classical problem 

rb 

%(•)] = / L(t,q(t),q(t))dt — ► min . 
•/a 

Theorem 3.1 (cf. [I'D If q is a minimizer of problem fljj, 
f/zen satisfies the fractional Euler-Lagrange equations: 

d 2 L(t,q, a D?q)+ t D%d 3 L{t,q, a D?q)=0. (5) 

The following definition is useful in order to introduce 
an appropriate concept of fractional conservation law. We 
recall that the classical Noetherian conservation laws are al- 
ways a sum of products (as assumed in (O) and that the frac- 
tional rule for differentiation of a product, in the sense of 
Riemann-Liouville, is enough complex (see e.g. [27]). With 
respect to this, our operator 9f (f, g) is useful. This operator 
was introduced in Hill and we refer the reader to this refer- 
ence for several illustrative examples and remarks. Here we 
just mention that the operator £) t a (f,g) has resemblances 
with the classical Poisson bracket (cf. Remark [34l l. 

Definition 3.1 (cf. |11|) Given two functions / and g of 
class C 1 in the interval [a,b], we define the following op- 
erator: 

®?{f,g)=faD?g-g t D%f, t€[a,b]. 
Remark 3.2 For a = 1, operator S> t a is reduced to 

@Hf,g)=faD}g-g t Dlf 

= fg+fg = jj(fg)- 

Remark 3.3 The linearity of the operators a Df and t D® im- 
ply the linearity of the operator 

Definition 3.2 (cf. Ifllll ) We say that Cf (t,q, a D"q) is ^frac- 
tional conservation law if and only if it is possible to write 
Cf in the form of a sum of products, 

C f (t,q,d) =£C? (t,q,d)-Cf(t,q,d) (6) 
i=l 



for some r 6 N, and for each i = 1 , . . . , r the pair C\ and Cf 
satisfy one of the following relations: 

9? (C? (t,q, a D?q) ,C? (t,q, a D?q)) = (7) 
or 

X {Cf (t,q, a D?q),Cj {t,q, a D?q)) = (8) 

along all the fractional Euler-Lagrange extremals (i.e. along 
all the solutions of the fractional Euler-Lagrange equations 
©). We then write 9 {C f (t,q, a D?q)} = 0. 

Remark 3.4 For a = 1 Ql and ([8]l coincide, and 

@{C{t,q, a D?q)} = 

is reduced to 

— {C(t,q(t),q(t))} =0^C(t,q(t),q(t)) =constant, 

which is the standard meaning of conservation law, i.e. a 
function C(t,q,q) preserved along all the Euler-Lagrange 
extremals q(t), t £ [a,b], of the problem. This implies that 
if (p(t),q(t)) is a solution to the classical Hamilton-Jacobi 
equations of motion, then C defines a conservation law of the 
Hamiltonian equations with Hamiltonian H if {H,C} = 
or {C,H} = 0, where {•,•} denotes the canonical Poisson 
bracket operator. In the more general fractional context, the 
Hamilton-Jacobi equations were recently derived in II18II191 . 

Definition 3.3 (cf. [11|) Functional © is said to be invari- 
ant under the one-parameter group of infinitesimal transfor- 
mations 

f = t + ex{t,q) + o{e), 

q(t)=q(t)+£^(t,q)+o(e), 
if, and only if, 

[ l ' L(t,q(t), ta D?q(t))dt 

= f L{t,q{t), T D?q{t))dt (10) 

Jt{ta) 

for any subinterval [?«,?/,] C [a,b]. 

Remark 3.5 Having in mind that condition dlOl) is to be sat- 
isfied for any subinterval [t a ,tj,] C [a,b], we can rid off the 
integral signs in ( 1101 ). This is done in the new Definition |4.3l 

The next theorem provides an extension of the classical 
Noether's theorem to Fractional Problems of the Calculus of 
Variations. 

Theorem 3.2 (cf. iflllP If functional (|4| is invariant under 
©, then 

[L(t,q, a D?q) ~ ad 3 L(t,q, a D?q) ■ a D?q]T{t,q) 

+ d 3 L{t,q, a D?q)-$(t,q) 

is a fractional conservation law (cf. Definition \3.2[ . 
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4 Main Results 



Using Theorem l3.2l we obtain here a Noether's Theorem for 
the fractional optimal control problems introduced in [ 3 ] : 



mm, 



(11) 



/[?(•),"(•)] = / L{t,q(t),u{t))dt 

J a 

a D?q(t) = <p(t,q(t),u(t)) , 

together with the initial condition q(a) = q a . In problem 
( flTi the Lagrangian L : [a,b] x W x R m -> M and the ve- 
locity vector (p : [a,/?] x M" x M m — > M" are assumed to be 
C 1 functions with respect to all the arguments. In agreement 
with the calculus of variations, we also assume that the ad- 
missible control functions take values on an open set of R™. 

Definition 4.1 Apair (#(•),«(•)) satisfying the fractional con- 
trol system a Dfq[t) = (p(t,q(t),u(t)) of problem dTTT) . t G 
[a, b], is called a process. 

Theorem 4.1 (cf. (13)-(15) of |3|) If (?(•),«(•)) « an o/?rf- 
»ifl/ process for problem (II lb . f/jen f/iere a co-vector 

function /?(•) smc/z f/jaf the following conditions hold: 

- the Hamiltonian system 

\ a D?q{t) = d 4 Jf(t,q(t),u(t),p{t)), 
\ t D«p(t) =dxJK , {t,q{f),u{t),p{t))\ 

- the stationary condition 
d 3 J^(t,q(t),u(t),p(t)) = 0; 

with the Hamiltonian J%? defined by 

J4? (t,q,u,p) =L(t,q,u)+p-(p(t,q,u) . (12) 

Remark 4.1 In classical mechanics, the Lagrange multiplier 
p is called the generalized momentum. In the language of 
optimal control, p is known as the adjoint variable. 

Definition 4.2 Any triplet (<?(•), «(•),/?(•)) satisfying the con- 
ditions of Theorem l4.1l wiil be called a fractional Pontryagin 
extremal. 

For the fractional problem of the calculus of variations 
(|4]l one has q>(t,q,u) = u =>• = L + p ■ u, and we obtain 
from Theorem [4j] that 

= u , 
t D%p = d2L, 

d 3 Jf = 0^p = -d 3 L =*> = - t D%d 3 L. 

Comparing the two expressions for iD^p, one arrives to the 
Euler-Lagrange differential equations q): d%L = — t D"d 3 L. 

We define the notion of invariance for problem dill) in 
terms of the Hamiltonian, by introducing the augmented func- 
tional as in [0: 

J[q(-),u(-),p(-)} 

= ( [Jf{t,q{t),u{t),p{t))-p{t)- a D t a q{t)]dt, (13) 



Remark 4.2 Theorem |4. 1 l is easily obtained applying the nec- 
essary optimality condition (0) to problem (113b . 

Definition 4.3 A fractional optimal control problem ( fTTT ) is 
said to be invariant under the e-parameter local group of 
transformations 



(14) 



t = t + ez(t,q(t),u{t), P (t))+o{e) , 
q(t)=q(t) + e%(t,q{t),u(t),p(t))+o{e), 
u(t) = u(t) + ea(t,q(t),u(t),p(t)) + o(e) , 
p(t)=p(t) + e£(t ) q(t),u(t),p(t)) + o(e) ) 

if, and only if, 

[Jf(t,q(t),u(t),p(t)) - p(t) ■ a D f a q(t)}dt 

= [^{t,q{t),u{t),p{t))~p{t)- a D t a q{t)]dt. (15) 

Theorem 4.2 (Fractional Noether's theorem) If the frac- 
tional optimal control problem dl lb is invariant under J 14b . 
then 

[JP - (1 - a) p[t) ■ a D t a q(t)] T - pit) ■ % (16) 

« a fractional conservation law, that is, 

@{[Jf - (1 - a) KO • aA a ?(01 t - p(0 • 4 } = 

a/ong aZZ the fractional Pontryagin extremals. 

Remark 4.3 For a = 1, the fractional optimal control prob- 
lem (fTTT l is reduced to the classical optimal control problem 

min, 



I[q(-),u(-)]= L(t,q(t),u(t))dt 

J a 

q(t)=(p(t,q(t),u{t)) , 

and we obtain from Theorem |4.2| the optimal control version 
of Noether's theorem [29|: invariance under a one-parameter 
group of transformations d 14b imply that 



C(t,q,u,p) =J?{t,q,u,p)T- p-t, 



(17) 



is constant along any Pontryagin extremal (one obtains d 1 7b 
from ( [TBI setting a = 1). 

Proof The fractional conservation law (116b is obtained ap- 
plying Theorem [32] to the augmented functional ( fT3l . □ 



Theorem 14.21 provides a new interesting insight for the 
fractional autonomous variational problems. Let us consider 
the autonomous fractional optimal control problem, i.e. the 
situation when the Lagrangian L and the fractional velocity 
vector (p do not depend explicitly on time t: 



min . 



(18) 



/[?(•),«(•)] = / L(q(t),u{t))dt - 

J a 

a D?q(t) = q>(q(t),u(t)) ■ 

Corollary 4.1 For the autonomous problem d 1 8b the follow- 
ing fractional conservation law holds: 



where is given by (fl2l) . 



3>{Jf-(l-a)p(t)- a D t a q(t)} = 0. 



(19) 
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Remark 4.4 In the classical framework of optimal control 
theory one has a = 1 and our operator Si coincides with 4 . 
We then get from ( fl9l ) the classical result: the Hamiltonian 
Jtf? is a preserved quantity along any Pontryagin extremal of 
the problem. 

Proof The Hamiltonian 3^ does not depend explicitly on 
time, and it is easy to check that dl8l) is invariant under time- 
translations: invariance condition ( 115b is satisfied with t = 
t + e, q(t) = q(t), u(t) = u{t) and p(t) — p(t). In fact, given 
that dt = dt, ( fT3T l holds trivially proving that s Df a q(t) = 
a D, a q{t): 

aDfq{t) 

1 / d \ " r'+ E 

= — — / (t+e-e) n - a - l q(e)de 

r(n-a) \dt) Ja+e V ' V ' 

= — — -( 4-V f {t-s) n ~ a - X q{s + e)ds 

r{n~a) \dt) J a 1 V ' 

= a D t a q(t + e)= a D t a q(t) 

= a D, a q{t). 

Using the notation in J 14b . one has T = 1 and ^ = a = £ = 0. 
Conclusion ([T9T l follows from Theorem |4.2| □ 



5 Illustrative Examples 

We begin by illustrating our results with two Lagrangians 
that do not depend explicitly on the time variable t. These 
two examples are borrowed from [3 §4.1] and [24, §3.1], 
where the authors write down the respective fractional Euler- 
Lagrange equations. Here, we use our Corollary 14. II to ob- 
tain new fractional conservation laws. 

Example 5. 1 We begin by considering a simple fractional 
problem of the calculus of variations (see 0] Example 1] 
and EU §3.1]): 

%(•)] = \ I (oD?q(t)) Z dt ^min, a>\. (20) 
2 Jo 2 

Equation (112b takes the form 

J^ = -lp 2 . (21) 
We conclude from Corollary [4j] that 

y(l-2a) (22) 
is a fractional conservation law. 



Example 5.2 Let us now consider the following fractional 
optimal control problem [3 §4.1]: 

I\q(,-)] = \j [q 2 {t)+u 2 {t)]dt ^min, (23) 
D?q(t) = -q{t)+u{t), 

under the initial condition q(0) = 1. The Hamiltonian Jff 
(112b has the form 

= ^(q 2 +u 2 )+p(-q + u). 
From Corollary 14- 1 l it follows that 

^(q 2 + u 2 )+ap{-q + u) (24) 

is a fractional conservation law. 

For a = 1 , the fractional conservation laws d22b and J24b 
give conservation of energy. 

Finally, we give an example of an optimal control prob- 
lem with three state variables and two controls (n = 3, m = 
2). The problem is inspired in [12, Example 2]. 

Example 5.3 We consider the following fractional optimal 
control problem: 

rb 

I (wi(f) 2 + u 2 (t) 2 ) dt — >min, (25) 

J a 

a D?q\{t) =wi(f)cos(g 3 (0), 

a D?q 2 {t) =a 1 (f)sinfe(f)), (26) 
aDfqiit) =u 2 (t). 

For a = 1 the control system d26b serves as model for the 
kinematics of a car and (I25b - d26l) reduces to Example 2 of 
H121 . From Corollary !4.1l one gets that 

u\ + ul+pi (miCos(<? 3 ) - (1 -a) a D?q\) 

+p 2 (mi sin(<? 3 ) - (1 - a) a D? q 2 )+P3 ("2 - (1 - a) a D^q 3 ) 

is a fractional conservation law. 

Main difficulty of our approach is related with the com- 
putation of the invariance transformations. To illustrate this 
issue, let us consider problem (1111) with 

L(t,q, u) — L(t, u) , <p(t,q, u) = (p(t, u) . 

In the classical case, since q does not appear both in L and (p, 
such a problem is trivially invariant under translations on the 
variable q, i.e. condition d 1 5b is verified for a = 1 with I = t, 
q(t) =q(t)+e, u(t) = u(t) and p(t) = p(t). In the fractional 
case this is not in general true: we have dt =dt, but condition 
( TT5T l is not satisfied since g D"q(t) = a D"q(t) + a D"e and the 
second term on the right-hand side is in general not equal to 
zero (Remark |22| . 
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6 Conclusions 

The fractional Euler-Lagrange equations are a subject of stron; 
current study lfT1l5l[8l l9l[T4lll5lll7ll25ll26l because of its nu- 
merous applications. In ifTTH a fractional Noether's theorem 
is proved. 

The fractional Hamiltonian perspective is a quite recent 
subject, being investigated in a serious of publications 
l^[TOll2T1l22llB1l24ll28l . One can say, however, that the frac- 
tional variational theory is still in its childhood. Much re- 
mains to be done. This is particularly true in the area of frac- 
tional optimal control where results are a rarity. The main 
study of fractional optimal control problems seems to be 0> 
where the Euler-Lagrange equations for fractional optimal 
control problems (Theorem 14. lb are obtained, using the tra- 
ditional approach of the Lagrange multiplier rule. Here we 
use the Lagrange multiplier technique to derive, from the 
results in 1111 . a new Noether-type theorem for fractional 
optimal control systems. Main result generalizes the results 
of 11291 . As an application, we have considered the fractional 
autonomous problem, proving that the Hamiltonian defines 
a conservation law only in the integer case a = I. 
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